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Abstract 

The  rendezvous  and  formation  problem  is  a  significant  part  for  the  unmanned  aerial  vehicle  (UAV)  autonomous 
aerial  refueling  (AAR)  technique.  It  can  be  divided  into  two  major  phases:  the  long-range  guidance  phase  and  the 
formation  phase.  In  this  paper,  an  iterative  computation  guidance  law  (ICGL)  is  proposed  to  compute  a  series  of 
state  variables  to  get  the  solution  of  a  control  variable  for  a  UAV  conducting  rendezvous  with  a  tanker  in  AAR. 
The  proposed  method  can  make  the  control  variable  converge  to  zero  when  the  tanker  and  the  UAV  receiver  come 
to  a  formation  flight  eventually.  For  the  long-range  guidance  phase,  the  ICGL  divides  it  into  two  sub-phases:  the 
correction  sub-phase  and  the  guidance  sub-phase.  The  two  sub-phases  share  the  same  iterative  process.  As  for  the 
formation  phase,  a  velocity  coordinates  system  and  a  nonlinear  decreasing  function  which  represents  control 
acceleration  are  designed  to  make  the  speed  of  the  UAV  consistent  with  that  of  the  tanker.  The  simulation  results 
demonstrate  that  the  proposed  ICGL  is  effective  and  robust  against  wind  disturbance. 

Keywords:  Autonomous  aerial  refueling;  Aerial  rendezvous;  Formation  control;  Guidance  law;  Unmanned  aerial 
vehicle 


1.  Introduction1 
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With  the  rapid  development  of  UAV  technologies,  UAVs  are  being  used  to  carry  out  various 
missions  like  high  altitude  surveillance  and  reconnaissance,  long  distance  military  strikes,  etc.,  some  of 
which  require  UAVs  continuously  stay  in  the  air  for  a  quite  long  time.  To  this  end,  autonomous  aerial 
refueling  (AAR)  becomes  a  key  issue  to  be  addressed  for  these  applications  of  UAVs.  There  are  four 
main  tasks  to  be  accomplished  sequentially  during  AAR  of  UAVs:  UAVs  rendezvous,  formation 
maintaining,  pipeline  docking,  and  refueling.  This  paper  mainly  focuses  on  presenting  a  guidance  law 
for  the  rendezvous  and  formation  maintaining  processes.  For  decades,  a  lot  of  achievements  have  been 
made  to  develop  guidance  technologies  for  the  UAV  rendezvous  problem  [1].  One  of  the  most  mature 
technologies  is  the  classical  proportional  navigation  guidance  (PNG)  [2].  After  that,  many  improved 
proportional  navigation  guidance  (IPNG)  schemes  have  been  evolved  based  on  PNG  [3,  4,  5].  To 
further  enhance  the  guidance  performance,  modern  control  theories,  such  as  sliding-mode  control  [6], 
differential  geometric  method  [7],  neural  network  [8],  and  so  on,  are  also  employed  for  the  design  of 
UAV  rendezvous  guidance  laws.  For  formation  maintaining,  there  are  several  existing  methods  such  as 
feedback  control  presented  in  Ref.  [9],  in  which  a  full-state  linearization  via  a  dynamic  feedback 
controller  is  designed  for  controlling  two  robots  in  a  leader-follower  configuration.  In  Ref.  [10],  a 
synchronized  position  tracking  controller  is  incorporated  in  formation  flight  control  for  multiple  flying 
wings.  Ref.  [11]  proposed  a  new  approach  of  hybrid  supervisory  control  for  the  leader-follower 
rmation  problem.  The  hybrid  supervisory  control  approach  provides  a  tractable  framework  for  hybrid 
esis  of  formation  control.  Within  this  framework,  a  new  method  of  abstraction  based  on  polar 
titioning  of  the  state  space  is  introduced.  Ref.  [12]  presented  an  iterative  guidance  method  for 
launch  vehicles.  In  this  method,  the  guidance  for  a  launch  vehicle  is  formulated  as  an  optimal  control 
problem,  in  which  the  transient  state  of  the  vehicle  is  taken  as  the  initial  value  and  the  target  point  as 
the  terminal  constraint.  The  objective  function  is  to  minimize  the  flight  time  of  the  vehicle  moving 
from  the  current  position  to  the  target  point.  During  the  whole  flight  process  of  the  vehicle,  for  each 
time  interval,  the  control  solution  and  the  corresponding  flying  trajectory  are  obtained  by  solving  the 
established  guidance  equations.  Through  a  repeated  iterative  computation,  the  launch  vehicle  is 
eventually  guided  to  the  target  point  and  satisfies  its  predefined  state.  This  method  has  been  studied  for 
the  guidance  and  control  design  of  launch  vehicles  and  ballistic  missiles  [13,  14,  15,  16].  However, 
these  approaches  involve  heavy  computation  load  and  are  difficult  for  practical  engineering 
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applications. 

In  this  paper,  a  novel  iterative  computation  guidance  law  (ICGL)  is  proposed  for  a  UAV  to  perform 
rendezvous  and  formation  maintaining  with  a  tanker  in  the  AAR  process.  In  the  ICGL,  the  rendezvous 
process  of  a  UAV  with  a  tanker  in  AAR  is  divided  into  two  major  phases:  the  long-range  guidance 
phase  and  the  formation  phase.  In  each  phase,  the  ICGL  computes  the  relative  state  parameters  between 
the  UAV  and  the  tanker  in  real-time  to  obtain  the  specific  control  variable.  Furthermore,  the  ICGL 
divides  the  long-range  guidance  phase  into  two  sub-phases  which  share  the  same  iterative  process:  the 
correction  sub-phase  and  the  guidance  sub-phase.  In  order  to  correct  the  actual  distance  error  between 
the  UAV  and  the  tanker,  an  estimated  distance  term  with  a  coefficient  is  added  to  the  model  in  the 
formation  phase.  By  adjusting  the  coefficient,  the  error  can  be  eliminated.  The  ICGL  proposed  in  this 
paper  is  totally  different  from  the  above  mentioned  iterative  guidance  method  based  on  the  optimal 
control  theory.  In  the  ICGL,  a  UAV  approaches  to  a  tanker  along  a  smooth  arc  trajectory  designed  by  a 
geometric  method.  Thus  the  iteration  is  based  on  the  relative  position  between  the  UAV  and  the  tanker 
in  each  time  interval.  The  designed  arc  trajectory  strategy  has  an  advantage  which  makes  the  control 
vector  perpendicular  to  the  velocity  of  the  UAV.  Therefore,  it  is  easy  to  be  realized  in  practical 
engineering  perspective. 

The  remainder  of  the  paper  is  organized  as  follows.  In  Section  2,  the  problem  of  rendezvous  for  AAR 
is  formulated  mathematically.  In  Section  3,  the  ICGL  algorithm  is  developed  for  a  UAV  to  perform 
rendezvous  and  formation  maintaining  with  a  tanker  in  AAR.  In  Section  4,  simulations  are  performed 
to  verify  the  effectiveness  of  the  proposed  ICGL  by  comparing  it  with  the  nonlinear  guidance  (NG) 
method.  Then  simulation  of  wind  disturbance  injection  is  conducted  to  demonstrate  the  robustness  of 
the  ICGL.  The  conclusion  remarks  are  given  in  Section  5. 


2.  Problem  description 


An  illustrative  diagram  for  the  rendezvous  process  of  a  UAV  with  a  tanker  for  aerial  refueling  is 
given  in  Fig.l.  The  research  objective  of  this  paper  is  to  guide  the  UAV  near  the  tanker  and  then  the 
UAV  keeps  a  formation  flight  with  the  tanker. 


r 


An  illustrative  rendezvous  process  of  a  UAV  with  a  tanker 

The  system  can  be  described  in  the  Earth-fixed  inertial  frame  Oc-RF  shown  in  Fig. 2  with  its  origin  at 
the  initial  location  of  U,  which  stands  for  the  UAV.  The  orientation  of  the  X-axis  is  consistent  with  the 
initial  velocity  of  the  UAV.  The  Z-axis  is  downward  and  perpendicular  to  the  X-axis  in  the  vertical 
plane.  The  Y-axis  is  determined  according  to  the  right-hand  rule.  The  transitional  state  of  U  during  the 
rendezvous  process  is  shown  in  Fig. 3. 


Fig.2  The  unparallel  situation  Fig. 3  The  state  of  the  UAV  during  the  guidance  process 

The  position  vector  of  U  is  denoted  as 
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The  velocity  vector  of  U  is 

The  motion  of  U  can  be  described  by  the  following  kinematics  equations  [16] 

fv,  ■  %  i7  1  CTK  T 
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u  is  the  angle  between  the  velocity  of  U  and  the  X-axis,  y  is  the  angle  between  the  guidance  surface 
m,  which  is  determined  by  the  guidance  acceleration  and  the  velocity  of  U,  and  the  plane  as 

shown  in  Fig. 3. 

In  Fig.2,  F  is  the  tracking  point  for  U  to  keep  a  formation  flight  with  the  tanker.  The  location  of  F 
can  be  set  according  to  the  current  position  of  the  tanker  with  a  deviation  value  Id  as  shown  in  Fig.2. 
The  position  vector  of  F  can  be  expressed  as 


»f  "ft 
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The  velocity  vector  of  F  is 


V-f 


where  is  the  position  vector  of  the  tanker  and  is  the  velocity  vector  for  the  tanker. 
The  motion  of  F  can  be  described  by  the  following  kinen 


3.  Guidance  law  design 
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With  the  conditions  discussed  above,  there  exist  two  situations  for  the  relationship  of  the  initial 
velocity  between  U  and  F:  unparallel  as  shown  in  Fig.2  and  parallel  as  shown  in  Fig. 4. 

Since  the  parallel  situation  i^an  exceptional  case  of  unparallel,  the  unparallel  situation  is  mainly 
discussed  in  this  paper.  The  rendezvous  process  can  be  divided  into  two  phases:  the  long-range 
guidance  phase  and  the  formation  phase. 

3.1  Long-range  guidance  phase 

The  long-range  guidance  phase  is  defined  from  the  initial  instant  to  the  instant  when  the  distance 
between  U  and  F  is  equal  to  or  smaller  than  a  given  value.  When  vu  and  Vf  are  not  parallel,  the 
long-range  guidance  phase  should  also  be  further  divided  into  two  sub-phases:  the  correction  sub-phase 
and  the  guidance  sub-phase. 


,u 


Correction  sub-phase 


The  correction  sub-phase  is  to  adjust  the  direction  of  to  parallel  with  the  velocity  of  the 
projective  point  P  within  the  Qm  ^  Yn  plane  firstly.  The  geometric  relationship  between  U  and  P  in  the 
Qf  Xi  Yf  plane  is  shown  in  Fig.5. 
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plane 


Fig.5  The  geometric  relationship  in  the  two-dime 


.ensional 


In  Fig.5,  r  is  the  distance  between  U  and  P;  a  is  the  angle  of  line-of-sight;  t?  is  the  advance  angle 
of  U\  a  is  the  lateral  acceleration;  ex  is  the  unit  vector  of  the  X-axis  and  er  is  the  unit  vector  of  the 
relative  distance  vector.  The  strategy  used  in  this  paper  to  guide  U  from  the  initial  position  to  the 
location  of  P  is  an  arc  orbital  approximation.  U  can  move  along  the  arc  trajectory  h  from  its  current 
position  to  the  target  location  by  applying  a  normal  corrective  overload  on  i 
The  guidance  angle  of  U  satisfies  the  constraint 


where  the  radius  of  h  can  be  calculated  by 


2  'ffri  r; 


where  tj  ■  Cf  -  ?,  in  which  q  m  arcvus  (,  V1!  ,  )• 

Wjc  1  Sr 

According  to  the  analysis  above,  the  control  acceleration  can  be  calculated  by 

The  iterative  process  throughout  the  correction  sub-phase  is  shown  below  in  detail. 
(1)  Initialization 

rfers  at  the  initial  instant 


(4) 


Initialize  the  followii 
The  position  vector 


The  velocii 
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ctor  of  U 


ing  parameters 
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The  position  vector  of  P 

The  velocity  vector  of  P 

The  guidance  angle  of  U 

The  relative  distance  vector  between  U  and  P 

The  magnitude  of  the  distance  between  U  and  P 
The  angle  of  line-of-sight 
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The  advance  angle  of  U 


The  magnitude  of  the  normal  correction  acceleration 


(2)  Iterative  computation 


In  order  to  get  a  series  of  solutions  for  the  control  variable  g  using  the  iterative  method,  the 
differential  equations  in  Eq.  (1)  and  Eq.  (3)  should  be  discretized  by  using  the  following  method 

- 3 - 

By  this  way,  the  guidance  process  can  be  divided  into  multiple  discrete  sub-processes  as  shown  in 
Fig.6.  Each  sub-process  is  regarded  as  a  straight  line  motion  with  time  interval 


Fig.  6  The  evolution  of  the  UAV  velocity 


From  the  discrete  sub-processes,  the  value  at  each  instant  can  be  calculated  step  by  step.  Let  i  (7=1,  2, 
3. . .)  represent  the  7-th  iterative  step;  the  iterative  computation  procedure  can  be  illustrated  as  follows. 
Step  1:  Set  i= 1  and  /  —  0; 

Step  2:  Calculate  the  position  vector  of  U  at  t  instant 


where 


Step  3:  Calculate  the  position  vector  of  P  at  instant 


wh 


Step  4:  Update  the  velocity  vector  of  P  at  q  instant 


Step  5:  Calculate  the  relative  distance  vector  between  U  and  P  at  frE  instant 


Vy  ($  " 


Step  6:  Calculate  the  magnitude  of  the  distance  between  U  and  P  at  i\  instant 


Step  7:  Calculate  the  angle  of  line-of-sight  c  at  tE  instant  with  Eq.  (18) 
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Step  8:  Update  the  guidance  angle  of  U  at  instant 

Step  9:  Calculate  the  velocity  vector  of  U  at  t  instant  by  Eq.  (1) 

where 


lJv:  'i ‘ 
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^■v^-fhcfU 

*W  ■  0 

Step  10:  Calculate  the  advance  angle  of  U  at  £,■  instant 

Step  1 1 :  Calculate  the  magnitude  of  the  normal  correction  acceleration  at  £,■  instant  using  Eq.  (4) 

vz  2v£ 


m 

Step  12:  Calculated  the  corrected  angle  of  the  velocity  vector  for  U  at  £,■  instant 


Step  13:  Calculated  the  angle  between  the  velocity  vectors  of 


Step  14:  If  y£6|J  is  smaller  than  a  given 
computation  and  go  into  the  guidance  sub-ph; 


?  U  and  P  ai 


W- 

value,  0.0005  rad  in  this  p 
ase;  otherwise,  let  I  ■  I 


■—■it 

'  1  "  'W 


it  £[  instant 

arms 


is  paper,  terminate  the  iterative 
*1,  and  go  to  step  2. 


3.1.2  Guidance  sub-phase 

At  the  end  of  the  correction  sub-phase,  the  position  relationship  between  U  and  F  as  well  as  P  in 
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Fig.7  The  position  relationship  after  the  correction  sub-phase  Fig. 8  The  new  relationship  in  Og-RF 
The  state  parameters  at  this  moment  can  be  denoted  as  follows: 

The  position  vector  of  U 

Pw  ■  (A; 

The  velocity  vector  of  U 

yvi  "  £%;# t  hasn’t 

The  position  vector  of  P 

Vvt  ■ 

The  velocity  vector  of  P 
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The  position  vector  of  F 
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Vft  ■(%%%) 


where  ■  v?i,.  fit-.-  ■  zr-  are  the  required  position  components. 

The  velocity  vector  of  F 

^  ■  V>f 

In  order  to  use  the  iterative  computation  algorithm  similar  to  that  used  in  the  correction  sub-phase,  a 
coordinate  system  Og-RF  is  created  with  its  origin  at  the  current  location  of  U.  The  orientation  of  the 

X-axis  is  consistent  with  the  current  velocity  direction  of  U.  The  Z-axis  is  downward  and  perpendicular 

to  the  X-axis  in  the  vertical  plane.  The  Y-axis  can  be  determined  according  to  the  right-hand  rule  as 
shown  in  Fig.7.  The  coordinate  transformation  from  Oc-RF  to  Og-RF  is  obtained  by 

G  ■%  ■  C  -h  £?  (5) 

where  C  is  a  3^1  column  vector  in  Oc-RF ,  R  is  the  rotation  matrix  from  Oc-RF  to  Og-RF  which 
can  be  calculated  as  follows: 

SQS  &  ffffo  F  0] 

CPS  P  Q|  (6) 

0  0 


and  D  is  a  3  ^  L  displacement  vector: 


ri  ■  1  jy n  ,■  ■- 


(7) 


Fig. 8  indicates  the  new  relationship  between  U  and  F  in  Og-RF.  For  the  guidance  sub-phase,  the 
initial  conditions  for  iterative  computation  at  ^  instantjjg^cpmputed  as  follows: 

The  position  vector  of  U 


J 


$ 


t  are  comp 

r 


The  velocity  vector  of  U 

The  position  vector  of  F 

The  velocity  vector  of  F 

The  angle  between  the  guidance  surface  tr.  and  the  OXY  plane 
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The  iterative  computation  process  is  the  same  as  that  in  the  correction  sub-phase. 

3.2  Formation  phase 

The  formation  phase  is  to  consider  how  to  keep  the  UAV  and  the  tanker  flying  in  a  formation  for  the 
preparation  of  the  docking  of  the  refueling  probe  of  the  UAV  with  the  drogue  of  the  tanker.  When  the 
distance  between  the  UAV  and  the  tanker  is  less  than  a  specified  value,  a  velocity  coordinate  system 
Ov-RF  is  created  as  shown  in  Fig. 9.  Ov-RF  is  used  to  measure  the  target  location  for  computing  the 
control  variable.  Fig.9  shows  the  position  relation  between  U  and  F  in  Ov-RF  during  the  formation 
phase.  In  order  to  guide  the  UAV  to  the  formation  point,  the  direction  of  the  control  acceleration  cf 
can  be  obtained  according  to  the  velocity  coordinate  system  as  shown  in  Fig.  10.  The  magnitude  of  at- 
is  determined  according  to  the  coordinate  value  of  F  in  Ov-RF  as 

■  A 

where  A  is  a  constant  and  i5  is  the  angle  between  the  UAV  velocity  and  the  X-axis  of  Ov-RF. 


(8) 
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Fig.9  The  UAV  position  relative  to  the  formation  point  Fig.  10  The  direction  of  the  control  of  the  UAV 

Besides,  the  formation  flight  phase  in  AAR  requires  the  speed  of  the  UAV  be  consistent  with  th 
the  tanker.  For  this  reason,  the  control  acceleration  along  with  the  direction  of  the  UAV 
should  be  utilized  to  control  the  speed  of  the  UAV.  The  magnitude  of  depends  on  ti 
distance  between  U  and  F  as  well  as  their  current  speeds.  In  order  to  compute  the  exact 
the  distance  for  U  to  move  is  given  by 

r* 


'  8  '3  'BTBffrl  (—  1  ■Hf'f 


where  the  term  of  "k  ■  r"  represents  an  estimated  distance,  in  which  "ft"  is  an  estimation  coefficient. 
According  to  the  Newton’s  theorem,  the  flying  time  for  distance  s  can  be  calculated  as 


o 


The  magnitude  of  the  control  acceleration  can  be  calculated  by 

In  summary,  the  rendezvous  process  can  be  described  by  the  flow  chart  shown  in  Fig.  11. 


V*  4-  Vf 


(10) 
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Fig.  1 1  The  flow  chart  for  the  ICGL  approach 


4.  Simulations  and  analysis 

In  order  to  verify  the  effectiveness  of  the  ICGL  for  the  rendezvous  and  formation  maintaining 
problem  in  AAR,  a  simulation  is  performed  here  using  the  proposed  ICGL  by  comparing  with  the  NG 
approach  presented  in  Ref.  [17]  to  demonstrate  the  effectiveness  of  the  ICGL. 

In  this  simulation,  assume  that  the  initial  conditions  of  the  UAV  and  the  tanker  are  as  follows: 
Position  vector  of  the  UAV 
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Pv'-k  3  ■  3 

(lOOm/fcOftt 


CO,  0. 0J 

Velocity  vector  of  the  UAV 
Position  vector  of  the  tanker 

Wk 3  ■  Ofy  jffB,  Eja)  ■  CS  0  0  0  m ,  -  S  0  0  0  ra ,  - 1 0  0  Qm3 

Velocity  vector  of  the  tanker 

ffCi; 3  ■  JVf*!  vTj  svK .)  ■  0  m  /  s,  -  4  0  m/p  ,?3 

During  the  formation  phase,  the  relative  position  of  the  UAV  to  the  tanker  is  set  as 

cV  ■  SSnfr  g':  ■  SSixir^  ■  3Qtr. 

Simulations  for  the  rendezvous  and  formation  maintaining  process  of  the  UAV  with  the  tanker  are 
performed  by  using  the  ICGL  and  NG  methods,  respectively.  Fig.  12  shows  the  trajectories  of  the  UAV 


A 


Fig.  12  Trajectories  through  the  rendezvous  process  using  ICGL  and  NG 

The  control  accelerations  for  the  correction  and  guidance  sub-phases  are  shown  in  Fig.  13  and  Fig. 
14,  respectively.  Using  £  ■  }.'•  E"GJ[18],  we  can  compute  the  energy  consumptions  ■  llJ9i3 

for  the  ICGL  and  ■  1£7Q  for  the  NG,  respectively.  Obviously,  the  UAV  consumes  less  energy  to 
accomplish  the  rendezvous  process  by  using  the  ICGL  than  by  using  the  NG. 


Fig.  13  Control  history  in  the  correction  sub-phase 
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t(s) 

Fig.  14  Control  history  in  the  guidance  sub-phase 

The  altitude  change  histories  of  the  UAV  during  the  rendezvous  process  are  shown  in  Fig.  15  and 
Fig.  16  for  the  ICGL  and  the  NG.  By  comparison,  it  can  be  seen  that,  during  the  correction  sub-phase,  it 
takes  52  seconds  for  the  UAV  to  complete  its  velocity  direction  adjustment  using  the  ICGL  and  34.5 
seconds  using  the  NG.  For  the  guidance  sub-phase,  it  takes  103.4  seconds  to  complete  this  process 
using  the  ICGL  compared  with  211.5  seconds  using  the  NG.  In  general,  it  takes  217  seconds  in  total  to 
accomplish  the  whole  rendezvous  process  using  the  ICGL  and  300  seconds  using  the  NG.  It  is  obvious 
that  the  proposed  ICGL  has  better  performance  than  the  NG  for  the  rendezvous  problem. 
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Fig.  15  The  UAV  and  tanker  flight  altitude 
change  histories  using  the  ICGL 


Fig.  16  The  UAV  and  tanker  flight  altitude( 
change  histories  using  the  NG 
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Fig.  17  and  Fig.  18  show  the  change  histories  of  the  distance  and  the  angle  between  the  UAV  and  the 
tanker  during  the  whole  rendezvous  process  using  the  ICGL,  respectively.  It  can  be  seen  from  Fig.  17 
and  Fig.  18  that  the  UAV  eventually  keeps  a  formation  flight  with  the  tanker.  This  indicates  the 
effectiveness  of  the  proposed  ICGL  method. 


Fig.  17  Distance  between  the  UAV  and  the  tanker 


Fig.  18  The  change  of  angle  for  the  velocity  vectors 

To  test  the  robustness  of  the  proposed  ICGL  to  wind  disturbance,  we  add  a  step  wind  disturbance 
along  the  X-axis  with  amplitude  of  gTi-  ■  £  h  at  the  simulation  time  of  300s  in  the  formation 
maintaining  phase  as  shown  in  Fig.  19.  The  simulation  result  for  the  whole  rendezvous  process  is  shown 
in  Fig.  20.  Fig. 21  and  Fig.22  show  the  deviations  of  the  UAV  to  the  desired  distance  with  the  tanker  in 
the  X  and  Y  directions  in  the  presence  of  wind  disturbance.  It  can  be  seen  that  the  steady  state  error  is 
about  2m  in  the  X  direction  and  less  than  3m  in  the  Y  direction.  Fig. 23  shows  the  change  of  the  control 
acceleration  of  the  UAV  under  wind  disturbance.  The  simulation  result  indicates  that  the  proposed 
ICGL  has  good  robustness  against  wind  disturbance. 
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Fig.  19  A 
disturbance 
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step  wind  disturbance 
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Fig. 20  Trajectories  of  the  rendezvous  process  against  wind 
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guidance  phase  and  the  formation  phase,  in  which  the  long-range  guidance  phase  is  further  divided  into 
two  sub-phases:  the  correction  sub-phase  and  the  guidance  sub-phase. 

(2)  The  ICGL  solves  the  guidance  problem  of  each  phase  and  sub-phase  in  a  two-dimensional  space 
with  the  iterative  method  to  obtain  the  control  acceleration  for  the  UAV. 

(3)  Simulation  results  demonstrate  that  the  proposed  ICGL  is  effective  for  the  UAV  to  deal  with  the 
rendezvous  and  formation  problem  in  AAR  and  is  robust  against  wind  disturbance. 
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